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Abstract 

Quantum interference is shown to deliver a means of regulating the diffraction 
pattern of a thermal atomic beam interacting with two standing wave electric 
fields. Parameters have been identified to enhance the diffraction probability 
of one momentum component over the others, with specific application to Rb 
atoms. 
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The realization of the effect of quantum interference in manipulating and controlling 
atomic and molecular processes has opened up a new area of active investigation in quan- 
tum dynamics, called coherent control (CC) method [1]. The approach relies on the cre- 
ation of a non- stationary superposition state comprising of two or more eigenstates of the 
atomic/molecular system. This is achieved by suitably designing an ultrashort laser pulse 
so that only a selective number of eigenstates are involved in the interaction between the 
atom and the laser pulse. Dynamics of this superposition state is entirely different from 
that of the eigenstate in the sense that all the superposed eigenstates contribute coherently 
towards the final dynamical outcome. The coherence among the eigenstates towards a par- 
ticular dynamical observable can be regulated by altering the population of the eigenstates 
as well as their relative phases. Although this has been demonstrated lately, in the control 
of molecular photo-dissociation [2], optical racemization [3], scattering [4], molecular beam 
deposition [5] etc., there are far fewer applications of this new innovation. In this letter, 
we report on the coherent manipulation of the diffraction of a thermal atomic beam by 
two standing light waves thus, providing a useful means towards building mirrors or beam 
splitters for atomic interferometers. Several works works have already been devoted to the 
diffraction process, [6-9] e.g., the Kapitza-Dirac diffraction of atoms [8], Bragg diffraction 
of atoms [9], etc. Classical coherence (due to the presence of two or more optical fields) and 
the effect of the electric field configuration in modifying the diffraction pattern has also been 
studied by some authors [7]. It is thus timely to introduce a different method based on the 
amalgamation of the classical and quantum coherence where the classical coherence (CC) is 
created from the phase relation between two electric field and the quantum coherence (QC) 
is created out of a pre-preparation of the atomic beam such that every atom in the beam lies 
in a superposition state of its eigenstates. The relative contribution of the CC and QC play 
a crucial role in regulating the atomic momenta transfer and hence, in the atomic diffraction 
pattern. Diffraction of an atomic beam from a light grating involves the transfer of mo- 
mentum from the light field to the atoms and this depends on the excitation paths. Thus, 
altering excitation paths, much the way it is done in the CC of molecular photo-dissociation 
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[2], would lead to the control in the atomic diffraction. Below we describe the theory and 
provide the computational results for Rb atom. 

A thermal beam of Rb atoms is pre-prepared in a way that every atom in the beam 
exists in a superposition state comprising of two of its coherently related eigenstates (atomic 
coherence), i.e., 

= c^e-^lfa > +c 2 {t)e- iE *l n \4>2 > (1) 

where |0j > are the eigenstates of the atomic Hamiltonian, of energy Ei. This preparation is 
achieved by passing the thermal beam through a ultra-short laser pulse decribed elsewhere 
[11]. The thermal nature of the beam is characterized by its x-directional Gaussian distri- 
bution function g(x) = {jxo 2 ^)~ x ^e ^1, where the full width at half maximum (FWHM), 
2a x \J — ln{\/2) of the Gaussian function directly relates to the temperature of the beam. 
This is followed by the passage of the beam through two coherently related standing light 
waves (optical coherence) whose amplitudes, frequencies and phases can be varied. That is 
the composit standing electric field E(x,t) interacting with the superposed atomic state is 
given by 

E(x,t)=J2E J (u J ,e j ,x,t)f J (z)k (2) 
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where Ej(uj,0j,x,t) = fj(x)e UOjt + c.c. with c.c. representing the complex conjugate 
of the terms preceeding it, 9 F = \6 1 — 9 2 \ is the relative phase of the two standing waves 
(SW), Ej°\ ujj and kj are the peak amplitude, frequency and wave vector of the j-th SW 
of wave length Xj. This field is considered polarized along the z-direction (see Fig.l for the 
general configuration of the proposed scenario). The functions fj(x) and fj(z) are the light 
beam profiles of the j-th electric field along the x and z directions respectively. We assume 
the atomic momentum along z is sufficiently larger than that along x so that light forces 
along z can be neglected and the longitudinal velocity v z remains constant. This effectively 
means that we can replace z in fj(z) by z = Zi n + v z t, where Zi n (set as zero) is the nozzle 
position where the beam ejects from. Taking fi(z) = /^(z) = exp(— (z — z ) 2 /°"z)> we replace 
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it by f{t) = exp(—(t — t ) 2 /r 2 ), where r = a z jv z and to = z /v z represent the interaction 
region and the time when the SWs have the maximum intensity. Thus fj(z) represents 
the adiabatic entry and exit for the atomic beam. The x-dependence fj(x) is taken as 
fj(x) = 2cos(kjX + 9j) which means that j-th electric field consists of two components which 
counterpropagate each other. 

The excitation now takes place from the superposition state to the higher excited states 
of the atom (see Fig.l for the atomic transitions driven by the SW fields). We assume 
that the interaction time is sufficiently small so that the transverse kinetic energy absorbed 
by the atoms during the interaction can be neglected (Raman-Nath approximation) and 
also that the freqencies of the light fields are well detuned from the atomic resonances so 
that spontaneous emission can be neglected. Thus we have an effective three-level system 
interacting with the field given by Eq.(2). Within rotating wave approximation (RWA) the 
Schrodinger equation for this system is 



a i 

v° 3 ; 



if {t) E?=i SlWfje"?"* if(t) ZU " 
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//(/) Vj JY ' if(t)XUitoif>e > M° 3 7 

where £)™' n =< <f> m \ t 1 \4> n > E^p /h is the Rabi frequency, A™' n = l/h(E n — E m ) — ujj 
is the detuning of the j-th SW field and H the atomic dipole vector. Integration of the 
above equation followed by the Fourier transformation (FT) delivers di(k), d 2 (k) and a 3 (k), 
where tilda refers to the FT of the respective time function. The functions di(k), d2(k) 
and d 3 (k) jointly describe the diffraction pattern of the atomic beam. Within second order 
perturbative theory we obtain d\{k) and d 2 (k) as 

d,{k) = g(k)ci + f E ^I?e^~g{k - sk 3 ) - \ £ [c^W^I™ (4) 

+c4 2 'X' 3) 'F ) + caflf ' 3 >n^/^][ e ^^^ (fc _ s{kj + kjl)) 
+e si(o j -e jl) ~ g{k _ s{kj _ kji))] 

d 2 (k) = g(k)c 2 + | E ^' 2) l[ jT e^g{k - sk 3 ) - \ £ [c 2 nf^I^ (5) 

j,s j,j',s 
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+c 2 nf> 3) nf 3) i¥' j,) + c.nf^Q^h^'^^gik - s{k 3 + k f )) 

+e si ^-W~g(k - sikj - k r ))] 

where the initial conditions (r = 0) are di(k) = Cig(k), a 2 (k) = c 2 g(k) and a 3 (k) = 0. If 
the intensity of the SW fields are relatively low and/or the detunings from the resonances 
are relatively large, the excitation from the initial pre-prepared state given by Eq.(l) would 
be constrained significantly and the contribution of as(k) would remain negligible. This 
is often referred to as adiabatic evolution of the initial state [10]. In such case, the light- 
atom interaction is equivalent to the propagation of a scalar atomic wave through an optical 
potential and the deflection of the atomic de Broglie waves by light fields is exactly analogous 
to that of the light by a dielectric medium [10]. The summation index s has two values viz., 
+ve and -ve, whereas l[ 3 \ I%' 3 '\ I^' 3 '\ and 7^' y) are given by 

/<.<> = r me -<<" it 

J — oo 

J —CO J —CO 

i™ = f dt'me-^' 3) f m^dt" 

J —CO J —CO 

J —oo J —oo 

and 

J —oo J —oo 

and * represents complex conjugate. The function g(k) is the Fourier transform of g(x). 
The scattered wave function of the atom is a superposition of the Gaussian modulated plane 
waves with momenta hk=0, ±hk±, ±hk 2 , ±h(ki — k 2 ), ±h(ki + k 2 ), ±2hk\ and ±2hk 2 . 
Thus, in the coherent control method the momentum transfered from the field to the atom 
has components different than when there was no coherence. The momentum h(ki + k 2 ) 
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(h(ki — k 2 )) corresponds to the absorption of a photon from the +ki component of the SW 
followed by induced emission in the — k 2 (k 2 ) component of the SW and is the result of the 
optical coherence. 

The output diffraction probability which can be measured experimentally, given by 
\di(k) + a 2 {k)\ 2 represents a comb of images of the split atomic velocities v = hk/m, where 
m is the mass of the atom. This output depends on the controlled parameters ci, c 2 , 
8m, Of, t, a x , e[°\ E 2 ^ and the superposed states, \<j>\ > and \<p 2 > with 9m being the 
atomic phase which is the relative phase between c\ and c 2 . The superposition is created 
between \n,li,mi > and \n, l 2 ,m 2 > Rydberg states employing a two photon preparatory 
stage (see ref. [11]), where n refers to the principal quantum number of the ground state of 
Rb, li and mj are the angular and azimuthal quantum numbers respectively. The Rydberg 
states were evaluated following the quantum defect theory. Thus, the process of excita- 
tion to the state |0 3 > occurs through three different paths viz., (1) the path that ends 
at |03 > from \n, l\,m\ >, the probability of which is proportional to |ci| 2 , (2) the path 
that ends at |0 3 > from \n,l 2 ,m 2 >, the probability of which is proportional to |c 2 | 2 and 
(3) the path that ends at \<f> 3 > through the interference between the paths (1) and (2), 
the probability of which is proportional to \ciC 2 \cos6 M , where 6m can be controlled during 
the preparation of the superposition state. We chose the frequency of the SW fields as 
ujj = Min{\E 2 — Exl/h, \E 3 — Ei\/h, \E 3 — E 2 \/h} + Aj which gives the detunings defined in 
Eq.(3) as A™'™ = l/h(E n - E m ) -Min{\E 2 - E^/h, {Es-E^/h, \ E 3 - E 2 \ /%} - Aj where m 
and n run from 1 to 3. In the calculation below for Rb atomic beam we chose Ai = 53051.6 
cmr 1 and A 2 = 477464.8 cmr 1 which ensure adiabatic evolution of the initial state. This 
means that the first and the second SW fields are detuned by Ai and A 2 respectively from 
the minimum of the atomic resonances. 

As an example, we examine the diffraction where the beam is peaked about zero along 
the transverse direction with a FWHM=5.5 jim which means the initial transverse spread 
(FWHM) of the atomic wave vector FWHM=0.504 //m -1 . This gives high resolution in the 
diffraction pattern of Rb beam as shown in Fig.2 for |ai(/c)| 2 (label a), |a 2 (/i:)| 2 (label b), 
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\a>i(k)\ 2 + \a2(k)\ 2 (label c) and \di(k) + a 2 (/c)| 2 (label d). Time integration is performed over 
a time period t = -5 x FWHM t to t = 5 x FWHM t , where FWHM t = 2r^-ln(l/2) 
is the full width at half maximum of the gaussian time profile centered at t — 0. In the 
computation we have taken FWHM t = 5/is which means that the superposed atomic beam 
interacts with the standing wave fields for a total duration of 50 /is [12], during which the 
fields become maximum at t=25 jis. In other words, for a fixed longitudinal velocity [13] 
v z = 500m/ s the field is maximum at distant of 12.5 mm from the nozzle. In Fig. 2 we 
observed the splitting of the beam into a total of 9 velocity components, of which the k=0 
is the most peaked. For the standing wave fields with Ai = 0.897/im {k-y = 7.004/im~ 1 ) 
and A 2 = 1.45/im (k 2 = 4.338/im" 1 ) the beam splits at k=0, 11.332, -11.332, 2.666, - 
2.666, 14.008, -14.008, 8.676 and -8.676 /zm -1 . We observed that the incorporation of more 
coherence (optical) through the introduction of more SW fields will give rise to the splitting 
of the atomic beam into a wide range of velocity components. Qualitative consideration of 
the pattern reveals that they are in general agreement with the theory [14]. That is, the 
atom absorbs photon from the +kj components followed by induced emission in the —ky 
components of the standing wave where j and f can be 1 and/or 2. 

We compute the relative diffraction probabilities of any two different momenta com- 
ponents Pk'k" = P(k')/P(k") where P(k') and P(k") are the total diffraction probabili- 
ties corresponding to the atomic momenta Tik' and hk" respectively. This P(k') is com- 
puted by filtering out all the momenta components in the diffraction profile excepting 
the one centered at k' and then integrating over the entire momentum coordinate, i.e., 
P(k') = f\ai(k) + a 2 (k)\ 2 F(k, k')dk, where F(k,k') is the filter chosen properly. Consid- 
eration of the equations for ay(k) and a 2 (k) (Eqs. (4) and (5)) show that the changes in 
the control parameters can strongly affect Pk'k") m other words, we can control the degree 
of splitting of one of the atomic velocity component over the others. The following figures 
(Fig. 3-5) show the dependence of Pk'k" on atomic phase 9m (Fig.3), |ci| 2 /|c 2 | 2 (Fig. 4) and the 
optical phase 9p (Fig. 5). All these figures show significant control of the diffraction probabil- 
ity of a given momentum component over the others by varying the controlled parameters. 

7 



For example, Fig. 3 (label c) show that the diffraction probability of k=0 component of the 
atomic wave vector is nearly 58 times that of the k = 2k 2 component at 6 M = 0.57 radian, 
i.e., -P(O) ~ 58P(2/c 2 ). This gets altered at 6m = 3.74 radian when the diffraction proba- 
bility of the k = 2k 2 component become nearly 11 times that of the k=0 component. Take 
another example, Fig. 4 (label j), where we see that the diffraction probability of k — 2k 2 is 
roughly 12 times that of k — 2k\ for 9m = —3.2 radian which totally reverses in favor of the 
k = 2k\ component at 9 M = —0.45 radian when P{2k\) « 11P(2A; 2 ). Thus the atomic phase 
causes a significant control on the diffraction probabilities of the atomic beam for fixed |ci| 2 
and |c 2 | 2 (Fig.3). 

Consider now the control over the Rb beam where initially |ci| 2 x 100 precentage of the 
atoms lie in the state |5, 0,0 > and the rest in the state |5,2,0 >. Results for such case 
are shown in Figs. 4 where we plot the logarithms of the relative diffraction probabilities 
of two different atomic wave vector components for varying | Ci | 2 / 1 C2 1 2 • Results show near 
monotonic decrease of P(0)/ P(h-k 2 ), P(0)/P(h+k 2 ), P{0)/P{2k 2 ) and P(0)/P(2k 1 ) with 
the increase in |ci| 2 /|c 2 | 2 for up to |ci| 2 /|c 2 | 2 ~ 3.5, after which they become nearly constant 
(Fig.4a). Whereas P(h - A; 2 )/P(A; 1 + k 2 ), P(h - k 2 )/P(2k 2 ), P(h - k 2 )/P(2k 1 ), P(h + 
k 2 )/P(2k 2 ), P(ki + k 2 )/P(2ki), P(2k 2 )/ P(2k 1 ) increase monotonically with the increase in 
|ci| 2 /|c 2 | 2 for up to |ci| 2 /|c 2 | 2 ~ 3.5 and then become constant (Fig. 4 b, c, d). Finally, Fig. 5 
depicts the control of the relative diffraction probabilities of two different atomic wave vectors 
for varying optical phase 0p, where we again observe several-fold increase (or decrease) of 
diffraction probabilities of a particular atomic wave vector over the other. 

The present work takes into account the experimental issues, e.g., (1) the experimental 
noise is taken care of by considering a gaussian thermal beam and (2) the entry and the exit 
of the atomic beam in its interaction with the SW fields is described by a gaussian envelope 
in the calculation. The parameters in the present calculation have been chosen carefully so 
as to commensurate with the Rb atom and the assumptions that the spontaneous emission 
is negligible and Raman-Nath regime is valid. These assumptions do not pose any obstacle 
in the real experimental scenario, they mere confine the intensity and the frequency of the 
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SW fields within certain values. 

In conclusion, we have shown that one can achieve significant control on the diffraction 
probability of a thermal atomic beam through the introduction of the atomic and optical 
coherences. Although we have introduced atomic coherence by two-state superposition sce- 
nario, a more general approach should be the preparation of the superposition state by an 
arbitrary numbers of atomic eigenstates accessible to the preparatory electric field. This 
would immediately introduce more parameters and hence more control over the diffraction 
pattern. 
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Figure Captions 

Figure 1: Schematic of proposed control scenario. 

Figure 2: Atomic diffraction pattern associated with the initial superposition states 
|0! >= |5,0,0 > and |0 2 >= |5,2,0 >. Here | d | 2 = 0.8, |c 2 | 2 = 0.2, eJ 0) = 1 x 10 4 V/m, 
E { 2 ] = 1 x 10 4 V/m, FWHM X = 5.5/im, FWHM t = 5.0/is, 6 F = 4.0 radian, 6 M =4.Q 
radian, Ai = 0.897/im, A 2 = 1.45/zm and t = 0. The lebels 1, 2, 3, 4, 5, 6, 7, 8 and 9 
correspond to the atomic momenta 0, h{ki — k 2 ), —h(ki — k 2 ), —2hk 2 , —Ti{k\ + k 2 ), —2hki, 
2hk 2 , h(ki + k 2 ) and 2hk\ respectively. The curves are normalized. 

Figure 3: Relative integrated diffraction probabilities Pk'k" corresponding to the atomic 
momenta (a) hk' = 0, hk" = h - k 2 ; (b) Hk' = 0, hk" = h + k 2 ; (c) Hk' = 0, hk" = 2k 2 ; (d) 
Hk' = 0, hk" = 2k ± ; (e) hk' = k x - k 2 , hk" = k x + k 2 ; (f) hk' = k x - k 2 , hk" = 2k 2 ; (g) hk' = 
fci - k 2 , hk" = 2h; (h) hk' = ki + k 2 , hk" = 2k 2 ; (i) hk' = h + k 2 , hk" = 2h; (j) hk' = 2k 2 , 
hk" = 2k\\ for the initial superposition states |0i >= |5, 0, > and |0 2 >= |5, 2, > plotted 
against 6m- Other parameters are as in Fig. 2. 

Figure 4: Logarithm of the relative integrated diffraction probabilities ln{Pwk") plotted 
against |ci| 2 /|c 2 | 2 for the initial superposition states |0i >= |5, 0,0 > and |0 2 >= |5,2,0 >. 
Here the lebels a, b, c, d, e, f, g, h, i, j signify the same as in Fig. 2. Other parameters are 
as in Fig. 2. 

Figure 5: Relative integrated diffraction probabilities Pyk" plotted against 9 F for the 
initial superposition states \<pi >= |5,0,0 > and |0 2 >= |5, 2,0 >. Here the lebels a, b, c, 
d, e, f, g, h, i, j signify the same as in Fig.2. Other parameters are as in Fig.2. 
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FIGURES 




FIG. 1. 
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FIG. 3. 
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FIG. 4. 
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FIG. 5. 



17 



